arXiv:1506.07762v2 [quant-ph] 12 Oct 2015 


Sagnac interferometry with coherent vortex superposition states in exciton-polariton 

condensates 

Frederick Ira Moxley III/ Jonathan P. Dowling/ Weizhong Dai/ and Tim Byrnes^’ 

^Heame Institute for Theoretical Physics, Department of Physics & Astronomy, 

Louisiana State University, Baton Rouge, Louisiana 70803-4001, USA 
^College of Engineering & Seience, Louisiana Teeh University, Ruston, LA 71272, USA 
^New York University, 1555 Century Ave, Pudong, Shanghai 200122, China 
^ NYU-ECNU Institute of Physies at NYU Shanghai, 

3663 Zhongshan Road North, Shanghai 200062, China 
^National Institute of Informaties, 2-1-2 Hitotsuhashi, Chiyoda-ku, Tokyo 101-8430, Japan 

(Dated: October 14, 2015) 

We investigate prospects of using counter-rotating vortex superposition states in non-equilibrium 
exciton-polariton Bose-Einstein condensates for the purposes of Sagnac interferometry. We first 
investigate the stability of vortex-antivortex superposition states, and show that they survive at 
steady-state in a variety of configurations. Counter-rotating vortex superpositions are of potential 
interest to gyroscope and seismometer applications for detecting rotations. Methods of improving 
the sensitivity are investigated by targeting high momentum states via metastable condensation, 
and the application of periodic lattices. The sensitivity of the polariton gyroscope is compared to 
its optical and atomic counterparts. Due to the large interferometer areas in optical systems and 
small de Broglie wavelengths for atomic BECs, the sensitivity per detected photon is found to be 
considerably less for the polariton gyroscope than with competing methods. However, polariton 
gyroscopes have an advantage over atomic BECs in a high signal-to-noise ratio, and have other 
practical advantages such as room-temperature operation and robust design. We estimate that the 
final sensitivities including signal-to-noise aspects are competitive with existing methods. 


I. INTRODUCTION 

Atomic Bose-Einstein condensation (BEC) was first 
achieved in 1995 [Ilia, which inspired a massive world¬ 
wide research effort into investigating its properties. To¬ 
day, atomic BECs are being considered for applications 
in a variety of fields, such as quantum simulation [3], 
quantum information mH], and quantum metrology 0. 
Recently, it has become apparent that BECs can be made 
in a variety of different systems. In particular, a BEC of 
quasi-particle excitations in semiconductor systems con¬ 
sisting of exciton-polaritons have been observed 13 m. 
Exciton-polariton BECs are formed in semiconductor mi¬ 
crocavities [9], or planar dielectric Fabry-Perot optical 
cavities m, for example, and can be described as half¬ 
matter, half-light excitations that are typically created 
using laser injection methods [niii2- Exciton-polariton 
BECs are a non-equilibrium phenomenon that can be 
achieved at room temperatures for suitable materials that 
support the polaritons m- In materials where polaritons 
can be operated at room temperature, interactions tend 
to be weaker, however more recently the nonlinearities 
have been observed to be significant El- Despite the 
non-equilibrium nature of the polariton BEC, it shares 
many of the properties of an atomic BEC, such as super¬ 
fluidity, stable vortex formation, and high spatial coher¬ 
ence. 

Although it is now commonplace to realize atomic 
BECs, typically they need to be cooled to very low tem¬ 
peratures in the region of ^ 100 nK. Even for high- 
temperature superconductors, which may be viewed as 
a type of charged condensation m, the highest criti¬ 


cal temperature attained has been 138 K [16]. On the 
other hand, exciton-polaritons (or polaritons for simplic¬ 
ity) typically have very high critical temperatures, due 
to their exceedingly light effective mass (i.e., ^ 10“^ 
times the bare electron mass). Consequently, this opens 
the door for new quantum technologies based on room- 
temperature superfluidity (i.e., “polaritonics”) m, and 
avenues for novel light sources |9|. 

One particular proposal based on BECs that has gath¬ 
ered recent interest is the matter wave analogue of the 
Sagnac interferometer, also called a “quantum gyro¬ 
scope.” An intuitive way to view these are as a neutral 
species version of a superconducting quantum interfer¬ 
ence device (SQUID) [18]. SQUIDs are well-known to be 
the most sensitive magnetometers, with capabilities of 
detecting small magnetic fields in the region of 5 x 10“^^ 
T m- Due to the close analogy of a magnetic field of 
a charged particle and rotation in the Schrodinger equa¬ 
tion, one may envision a device based on BECs that is 
highly sensitive to rotations, instead of magnetic fields. 
For this reason, ultracold atomic gases have been con¬ 
sidered to be an outstanding candidate physical system 
for realizing such a Sagnac interferometer [inillT]. One 
advantage that atomic systems possess is the exquisite 
experimental control that is achievable, particularly for 
creating vortices in BECs [22]. Moreover, these ultracold 
atomic gases are highly coherent [23] , meaning that they 
are particularly suited for interferometry. However, from 
a practical perspective, one of the drawbacks of technolo¬ 
gies that utilize ultracold atomic gases is that in order to 
achieve BEC, experimental conditions need to be cooled 
to nanokelvin temperatures. This can prove quite expen- 
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sive, and difficult to achieve in a compact device design. 
In order to overcome this challenge, a natural extension of 
the idea is to use exciton-polariton BECs [24l [25] , which 
may condense at room temperatures EillTj. This of¬ 
fers a potential new configuration for achieving a novel 
Sagnac interferometer [28|. In this paper we investigate 
in detail the possibility of such a “polariton quantum 
gyroscope.” The type of novel measurement configura¬ 
tion investigated herein has the potential for providing 
a highly compact and stable Sagnac interferometer [29] . 
which can also be used as a seismometer, for example. 
Future applications of quantum gyroscopes include the 
detection of earthquakes m, gravitational waves [311132], 
sidereal time [33], inertial sensing [34] capabilities, and 
even the rotational rate of our universe [35] . 

Such an exciton-polariton gyroscope would be a hy¬ 
brid optical-matter wave Sagnac interferometer [36], as 
polaritons themselves consist of both photons and exci- 
tons. In our scheme, one of the necessary ingredients 
for this interferometer will be the optical excitation of 
vortices m in the polariton BEC. The direct excita¬ 
tion of vortices by optical means is a convenient way of 
producing the required vortex superposition state in the 
condensate, as techniques to generate the optical coun¬ 
terparts are well-established. Here, it should be pointed 
out that the counter-rotating vortex superpositon state 
differs from the vortex-anti vortex pair such as those ob¬ 
served in Ref. [38]. The investigation of various vor¬ 
tex states in polariton BECs has continued to gain much 
interest, with several works examining disorder-induced 
vortices [39] , vortex-antivortex pairs EH 133 SO], and vor¬ 
tex ring creation via rotation of Gaussian laser beams 
El, to name a few. One of the most efficient methods 
of vortex creation is using Laguerre-Gauss optical states 
which carry a well-defined external orbital angular mo¬ 
mentum (0AM) [42l[44] . However, for Sagnac interfer¬ 
ometry, a superposition of a vortex and an antivortex is 
required, which has not been achieved in polariton BECs 
to date to our knowledge. While the vortex-antivortex 
superposition has been achieved in atomic systems [45] . 
this remains an outstanding problem for polaritons, and 
is interesting from both a fundamental and applications 
perspective. 

The purpose of this article is to investigate the 
prospects of angular momentum superposition states 
in exciton-polariton BECs for potential applications in 
Sagnac interferometry. In Sec. [TT]we present and analyze 
a scheme for generating vortex-antivortex superpositions 
via coherent excitation with 0AM states of light. In the 
scheme, the polariton BEC is initially seeded with the de¬ 
sired 0AM superposition, then allowed to evolve freely 
via the open-dissipative Gross-Pitaevskii equation. This 
allows for the polaritons to approach their own dynami¬ 
cal equilibrium state, independent of the initial state of 
light. We then study methods to improve the sensitivities 
of the polariton BEC gyroscope in Sec. |ln| This includes 
using ring geometries to target high angular momentum 
states, and ways of generating this using metastable con- 



FIG. 1: An experimental configuration for producing vortex- 
antivortex superpositions in exciton-polariton Bose-Einstein 
condensates. A laser generates a TEMqo Gaussian beam pro¬ 
file. This light is passed through a spiral phase plate (SPP), 
generating an orbital angular momentum (0AM) state of 
light. Once the 0AM state of light has traveled through the 
Mach-Zehnder interferometer with a Dove prism (DP) and 
mirrors (M), an 0AM superposition state of light is made. 
The resulting 0AM superposition state of light can be used 
to generate an interference pattern in a polariton BEC matter 
wave within the distributed Bragg reflector (DBR) microcav¬ 
ity. This interference pattern is imaged by a detector such 
as a charge-coupled device (CCD), and can then be used to 
determine the Sagnac phase and therefore the rotational rate 
of the polariton BEC. 

densation. A second enhancement scheme is discussed, 
where periodic lattices are imposed in order to intro¬ 
duce high momentum components to the wavefunction. 
In Sec. El we discuss the sensitivities of the Sagnac 
phase for the vortex superposition state, and compare 
with other existing approaches. In particular, we will de¬ 
rive the Sagnac phase for the vortex superposition and 
ring geometry cases, and compare the performance to op¬ 
tical and atomic Sagnac interferometers. Finally, we will 
summarize our results and conclude in Sec. V. 


II. VORTEX-ANTIVORTEX SUPERPOSITIONS 
IN POLARITON BECS 

A. Initialization by Laguerre-Gauss modes 

We first describe how the vortex-antivortex superpo¬ 
sition is created in the exciton-polariton BEC. The ex¬ 
perimental setup for the creation of a vortex-antivortex 
polariton BEC is shown in Fig. [^ which is similar to the 
traditional optical (i.e., Mach-Zehnder) Sagnac interfer¬ 
ometer. Unlike the traditional Mach-Zehnder scheme, 
our polariton BEC interferometer is based on a princi¬ 
ple that employs 0AM superposition states of light. Our 
starting point is a source of 0AM states of light, which 
can readily be produced by established methods, for ex- 
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ample, by using spiral phase plates (SPPs). This may 
be used to create 0AM superpositions, which coherently 
excite the polaritons in a DBR structure (e.g., semicon¬ 
ductor). This produces the initial seed polariton conden¬ 
sate, from which the polariton BEG is reinforced by the 
standard process of exciton-polariton cooling and scat¬ 
tering. 

Let us denote the bosonic annihilation operator for an 
0AM mode with angular momentum number I and radial 
node number p as ai^p (see Appendix [A| for more details). 
Then the superposition state for the mode that we re¬ 
quire is aai^p + l3a-i^p, where \a\‘^ + |/3p = 1. This 0AM 
superposition state of light has two counter-rotating com¬ 
ponents which will induce an interference pattern in the 
polariton BEG. Once the characteristic interference of 
the 0AM superposition (consisting of 21 lobes) has been 
transfered to the polariton BEG, the Sagnac phase can 
then be determined from the angular velocity of the BEG 
cloud. We can obtain an ordinary vortex state here sim¬ 
ply by setting either a or f3 = 0 [46]. 

Many lasers emit light beams which approximate a 
Gaussian profile with I = p = 0. Using the TEMqo fun¬ 
damental transverse mode of the optical resonator and 
channeling it through a SPP, for example, can be used 
to generate 0AM states of light (see Appendix 0- Al¬ 
though these states are generally made with holographs, 
SPPs offer a simple and convenient way to obtain a pure 
0AM state ai^p with arbitrary integer 1. Moreover, once 
the pure state ai^p has been obtained, we can generate an 
0AM superposition state. One way to accomplish this is 
by using a Mach-Zehnder type of configuration as seen 
in Eig. where the Dove prism is used to transform the 
ai^p state into a a-i^p state to obtain {ai^p + a-i^p)lV^ 
at the real output of the Mach-Zehnder interferometer. 
If desired, one may choose the first beam splitter as an 
|a|^:|/3|^ beam splitter, and the second one as a 50:50 
beam splitter to obtain an arbitrary superposition of the 
type aai^p -h pa^i^p. 

Let us now consider the specific case as shown in Eig. 
where the two inputs of the first beam splitter are 
prepared in a coherent state in the LG mode ai^p and the 
vacuum state respectively. Hence, 

UlN = CLi^p. (1) 

Assuming that the first beam splitter being used is a 
symmetric 50:50 beam splitter, the 0AM state of light 
has transformed like 


UlN = CLi^p ^ — {hi^p -h Ci^p) 


( 2 ) 


where bi^p^ ci^p denotes the LG modes in the two arms of 
the interferometer. The phase shift in the upper arm of 
the Mach-Zehnder interferometer will cause the coherent 
amplitude to replaced by ai^p e^^ai^p^ and the Dove 
prism in the lower arm of the Mach-Zehnder interferom¬ 
eter will cause the coherent amplitude to be replaced by 
bi,p ^ Hence, the 0AM state of light seen in Eq. 


© becomes 

(bi,p + Q,p) ^ (e ^bi^p + C-i,p^ • (3) 

Einally, neglecting the imaginary output of the interfer¬ 
ometer, the second beam splitter reverses the transforma¬ 
tion and produces the 0AM superposition state of light, 
written as 


^iouT = ^ {cLi,p + , (4) 

which is our desired output. The superposition of ai^p 
and a-i^p states in Eq. 0 can be detected using many 
different techniques, which typically includes imaging the 
interference pattern of the light beams which constitute 
the 0AM superposition state. In particular once the po¬ 
lariton BEG is induced, one may use standard spatial 
and momentum space imaging methods to analyze the 
coherent condensate state in the polariton microcavity 
0 . 


B. Exciton-polariton BEG dynamics 

The laser injection methods of the previous section 
can be used to resonantly excite the polariton BEG. The 
phase of the light is directly imprinted on the polari¬ 
ton BEG system, hence the injected polaritons simply 
follow the same phase relation as the injected light. In 
this configuration alone, there is no difference of using 
the 0AM superposition states of light to the polariton 
BEG, as the polaritons inherit their vortex nature from 
the light. However, once the initial coherent excitation 
is induced, the 0AM superposition is turned off, and 
the microcavity system is excited by conventional polari¬ 
ton pumping. This can be performed by an off-resonant 
pumping where the excitons have a much higher energy 
than the polaritons, or by pumping at a high in-plane 
momenta (see Eigure lb in Ref. 0)- Either method cre¬ 
ates a large population of uncondensed polaritons from 
which the polariton BEG is replenished. It is important 
to note that apart form the initial seed polariton BEG 
which is induced coherently, the remaining population 
thereafter is not pumped with the 0AM superposition 
state. Hence in the second phase where the conventional 
polariton BEG pumping is used, any 0AM superposition 
state is a genuine steady state of the polariton BEG, and 
possesses its own dynamics. 

While both vortices and vortex-antivortex pairs have 
been shown theoretically and experimentally to be stable 
configurations of polariton BEGs [ 31113 , whether or not 
a counter-rotating vortex superposition is a stable config¬ 
uration is a non-trivial problem. Such vortex-antivortex 
superpositions have a more complex phase relation than 
standard vortex configurations, and have the possibility 
to spontaneously separate spatially, or relax to a zero mo¬ 
mentum state. To investigate this, we have studied the 
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vortex-antivortex superposition in a variety of different 
circumstances to evaluate its stability. 

The spatiotemporal evolution of the non-equilibrium 
exciton-polariton system is described by the open- 
dissipative Gross-Pitaevskii (dGP) equation, expressed 
as p%l l48] 


.d^p{x,t) 
* dt 


- +Vext{x) + g\ip{x,t)f 


where ' 0 (x, t) is a complex valued condensate order pa¬ 
rameter, g is the polariton-polariton interaction constant, 
m is the polariton mass, Vext{'^) is a spatially dependent 
trapping potential energy (see Figure 5 in Ref. [ 8 ]), P(x) 
is a spatially dependent pumping rate, 7 is the polariton 
loss rate, and 7 is the gain saturation iiniiin]. Through¬ 
out this paper we will assume that the time scale is in 
units of 2 ma/fi, where a is the length scale set by the 
experiment (e.g., size of trapping potential or pump spot 
diameter). 

We use the Generalized Finite-Difference Time- 
Domain (G-FDTD) scheme to solve the dGP equation. 
G-FDTD is an explicit method that permits an accurate 
solution with simple computation for solving Eq. ^ in 
multiple dimensions. To this end, the function '0(x, t) is 
first split into real and imaginary components resulting 
in two coupled equations. The real and imaginary com¬ 
ponents are then approximated using higher-order Taylor 
series expansions in time, where the derivatives in time 
are then substituted into the derivatives in space via the 
coupled equations. Finally, the derivatives in space are 
approximated using higher-order finite difference meth¬ 
ods. The G-FDTD has been successfully applied for 
solving both linear and nonlinear Schrodinger equations 
[5ll|52]. 


C. Stability of vortex-antivortex superpositions 


We now examine the stability of the vortex-antivortex 
superposition under the dynamical evolution of the dGP 
equation. In Fig. we show the spatiotemporal evolu¬ 
tion of the non-equilibrium polariton BEG vortex super¬ 
position state in a flat potential landscape Vext = 0- The 
initial condition was chosen as 

ip{x,y,t = 0 ) = '^ {uf=i,p=o + 7 =-i,p=o)- ( 6 ) 

where are the spatial Laguerre-Gauss wavefunctions 

defined in Appendix 0. The factor of ^ approx¬ 

imates the equilibrium density which the dGP equation 
converges to. Eor densities near equilibrium 




P{x,y) - 7 
V 


( 7 ) 


the imaginary terms on the right hand side of Eq. 
vanish and we may obtain a zeroth order estimate for 
solutions [ 53 ]. Starting from this approximation, the 
numerics should converge towards the true solution for 
the dGP equation, which will in general be different 
from the standard (conservative) Gross-Pitaevskii solu¬ 
tion. Throughout this study, we assumed that the wave- 
function in the ^-direction has an insignificant effect on 
the vortex dynamics in the x-y plane. In the limit 
P = 7 = 7 = 0, the condensate is conservative (i.e., 
the particle number is fixed and there is nothing enter¬ 
ing or leaving the condensate). With P — 7, 7 > 0, the 
condensate has both leakage and gain, but stability is 
maintained such that there is no net loss. 

Erom the initial condition given by Eq. (§, the sys¬ 
tem is time-evolved under Eq. (§ until a steady state 
is reached. Eig. illustrates the condensate density 
distribution |'0(x,^,t)p of the two-dimensional vortex 
non-equilibrium condensate in a flat potential at various 
times, and also illustrates the phase (i.e., the argument of 
'0(x, t) ) of the two-dimensional vortex non-equilibrium 
condensate in a flat potential at various times. Upon ex¬ 
amination of Eig. we can see that our counter-rotating 
vortex-antivortex superposition state is coherent and sta¬ 
ble as it is generated in the distributed Bragg reflector 
(DBR) microcavity outlined in Eig. Here, it should 
be pointed out that the vortex-antivortex superposition 
has the characteristic phase dependence of the lobes be¬ 
ing out of phase with each other by a factor of tt. The 
system quickly reaches steady state and little change is 
seen beyond t = 10. The stability of the counter-rotating 
vortex-antivortex superposition can be attributed to the 
same reason as for atomic BECs: both the vortex and an¬ 
tivortex require a density defect due to the phase velocity 
diverging at the vortex core. Thus it is energetically fa¬ 
vorable for both to share the same vortex core, effectively 
pinning them in the same location. 

In Eig. [^we also simulate the vortex-antivortex super¬ 
position for experimental parameters as given in Refs. 
[H i5i Eg. We see that generally the results are un¬ 
changed and the interference lobes are stable under dif¬ 
ferent parameters. As the vortex superposition state con¬ 
sists of lobes, the Sagnac phase can be determined from 
the rotation of the lobes about an axis of rotation. This 
will be discussed in more detail in Sec. lEl 

While Eig. [^suggests that the vortex-ant ivort ex super¬ 
position is stable, a realistic semiconductor system does 
not have a perfectly smooth potential Vext = 0. There 
are typically unavoidable sources of disorder, which may 
potentially affect the stability of the superposition states. 
To investigate this, we study the evolution of the vortex- 
antivortex superposition state in a disordered potential 
landscape. Upon examination of Eig. we can see that 
in the presence of a disordered potential, our vortex su¬ 
perposition state is coherent and stable as it is generated 
in the distributed Bragg reflector (DBR) microcavity out¬ 
lined in Eig. 0 Fig. i also illustrates the phase of the 
two-dimensional vortex non-equilibrium condensate in a 
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density phase 




FIG. 2: Time evolution of a polariton BEG initially seeded 
with a vortex-antivortex superposition for a flat potential 
I4cct(x) = 0. Timeframes for the (a) initial state t = 0; (b) in¬ 
termediate state t = 5; and (c) steady-state t = 10 are shown. 
Parameters used were P{x,y) = , where Pq = 2, 

ro = 5.35 , and g = ^ = f] = 1. Timescales are in units of 
2mo? I % where a is the lengthscale unit. 


disordered potential, which shows in general the same 
behavior as the flat potential case. We have tried this 
for several disorder types and strengths, and find that 
the vortex-antivortex superposition is stable in all cases. 
Again, upon examination of Fig. we can see that with 
experimental parameters taken from Refs. m [Ml [55] 
the results are largely unchanged. 


III. SENSITIVITY ENHANCEMENTS FOR 
POLARITON GYROSCOPY 

The results of Sec. [TT| suggest that counter-rotating 
vortex-antivortex superpositions can be coherently ex¬ 
cited and are stable in polariton BECs. In the context 
of gyroscopy, rotations are detected by changes in the 
lobe interference pattern as found in Figs. i-i Thus 
for sensitive detections one would like the lobes to be as 
sharply defined as possible, to be able to discriminate be¬ 
tween the original and rotated interference patterns. For 


density phase 



FIG. 3: Time evolution of a polariton BEG initially seeded 
with a vortex-antivortex superposition for a flat potential 
Vecct(x) = 0. Timeframes for the (a) initial state t = 0; (b) 
intermediate state t = 5; and (c) steady-state t = 10 are 
shown. Timescales are in units of /i/meV. The following pa¬ 
rameters were used in the simulations: Tig — 0.05, /ly = 1.0, 
and hr] = 0.1, where energies are in meV. 


vortex-antivortex superpositions in polaritons, this can 
be achieved using large angular momentum states I (see 
Fig. 2 of Ref. [56]). While this is possible in principle, 
there are practical limits on the high angular momentum 
Laguerre-Gauss modes that can be excited. Hence alter¬ 
native strategies for sensitivity boosts are desirable for 
the polariton gyroscope. 


A. Metastable condensation in ring geometries 

In general, higher sensitivities are achieved by the use 
of high momentum states due to their short wavelength, 
allowing for high resolutions. A variant of the vortex- 
antivortex superposition state is to use ring geometries, 
where the polariton condensate is placed in a circular 
narrow channel potential. In a BEG, high angular mo¬ 
mentum states are energetically disfavored, and typically 
multi-quantized vortices break up into many / = ±1 vor¬ 
tices [53] . This can be understood to result from the high 
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density phase 


density phase 



FIG. 4: Time evolution of a polariton BEG initially seeded 
with a vortex-antivortex superposition for a disordered po¬ 
tential. Timeframes for the (a) initial state t = 0; (b) inter¬ 
mediate state t = 5; and (c) steady-state t = 10 are shown. 
Parameters used were P{x,y) = , where Pq = 2, 

ro = 5.35 , and g = ^ = f] = 1. Timescales are in units of 
2mo? I % where a is the lengthscale unit. 


FIG. 5: Time evolution of a polariton BEG initially seeded 
with a vortex-antivortex superposition for a disordered po¬ 
tential. Timeframes for the (a) initial state t = 0; (b) inter¬ 
mediate state t = 5; and (c) steady-state t = 10 are shown. 
Timescales are in units of TilmeV. The following parame¬ 
ters were used in the simulations: Tig = 0.05, /ly = 1.0, and 
hr] = 0.1, where energies are in meV. 


momentum states that are near the vortex core which 
have a rapid phase variation. By eliminating the vortex 
core with the use of a ring potential, high angular mo¬ 
mentum states become ultra-stable, which is beneficial 
from a sensitivity point of view. 

In Fig. [^we show the stability of counter-propagating 
polariton superfluid currents in ring geometries. The po¬ 
tential that is used is a “mexican hat” form 

Ve.tir) = Vo(^-2^), ( 8 ) 

where r = Expanding around the minima of 

the channel, one obtains an approximate wavefunction 


= exp 




{r - rminf 


(9) 


where we have added a phase variation around the ring 
with angular momentum 1. We then start in the initial 


state 


tp{r, (f)) 


(/)) + tp-iir, (j)) 

V2 


( 10 ) 


and evolve the system using Eq. ([^ . Eigure shows our 
results for evolving I = ±1 and I = ±5 superpositions. 
We see that steady state configurations are reached with 
the superpositions being maintained for all times, in the 
same way as for the vortex-antivortex superpositions. 

While the above shows that with a suitable initializa¬ 
tion of the the polariton BEG with counterpropagating 
currents is a stable configuration, a question is how in 
practice such a state would be induced. The method of 
coherent excitations as discussed in Sec. imi is viable 
for low angular momenta superpositions, but for higher 
angular momenta superpositions this can be problematic 
[56] . An alternative method of forming the superposi¬ 
tions states is via metastable condensation. In Ref. m 
it was experimentally found that in a one-dimensional 
periodic lattice, polariton condensates at high momen¬ 
tum states could be generated. This was attributed to 
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convenience that 2 p{x,y) = + L) although this is 

not important for the purposes of our illustration, and 
other boundary conditions could be used just as effi¬ 
ciently. 

Our target configuration in this case is that we have 
two counterpropagating momentum modes ±ko. In a 
similar way as aforementioned, in the ideal case where 
only these two modes are initialized, this is a stable con¬ 
figuration of the polariton BEG. To examine what hap¬ 
pens with an imperfect initialization procedure, let us 
consider that we have in addition a component with zero 
momentum 

y,t = 0)= a=o + ^0 • (12) 

We compare two cases where the pump is spatially uni¬ 
form 


P{x,y)=Po (13) 



FIG. 6: Counterpropagating currents in ring geometries. 
Cases shown are (a) I = ±1, t = 0; (b) / = ±1, t = 10; 
(c) I = ±5, t = 0; (b) / = ±5, t = 10. Parameters used are 
7 = ?7 = Vb = 5'=l, Tmin = 5. A spatially dependent pump 
P{x,y) = Po|'0z(r, 0)1^, with Po = 2 is used. Timescales are 
in units of 2ma^ jh where a is the lengthscale unit. 


the tradeoff between relaxation from high energy states 
and the finite decay time of the polaritons. Here we 
show that it is possible to target particular momentum 
state superpositions using the metastable condensation 
dynamics and spatially dependent pumping profile. 

For large rings we may model the ring geometry as 
discussed above by a toroidal geometry where the x di¬ 
rection is the direction along the channel, and y is in the 
radial direction. Then the ring geometry imposes that 
the wavefunction is periodic in x 

'tp{x,y) ='tp{x + L,y) (11) 

where L is the circumference of the ring. We also take for 


and a spatially varying pump which has the same pump¬ 
ing profile as the target wavefunction density, which in 
this case is cos^ k^x. 

P(x, y) = PoT] cos^ kox + 7 . (14) 

In the case of Ref. EH, the spatially periodic pump mod¬ 
ulation occurs because of the metallic deposits on the sur¬ 
face of the DBR which physically blocks the laser pump. 
These metal deposits also have the effect of producing a 
periodic potential of the form 

Vextix, y) = Vo cos(2kox), (15) 


induced by the change in the photonic boundary condi¬ 
tions at the surface [58]. The factor of 2 arises because 
for TT state condensation as observed in Ref. [57j, the 
phase in the wavefunction changes by tt between adja¬ 
cent lattices. 

Fig.Q shows our results. We start in a state with an 


equal distribution of momenta, ^k=o = C/c=/co = ^k=-ko = 
l/v^S. Under the uniform pumping scheme Eq. (13), the 
steady state evolves towards dominantly the zero mo¬ 
mentum mode. There is a spatial dependence to the real 
space wavefunction according to the potential Eq. (15), 
but the phase is the same across all the sites, indicat¬ 
ing that it is dominantly a zero momentum mode. This 
is explicitly seen by the Eourier space decomposition in 
Fig. 0b). Eor the periodic pumping Eq. ( [l4| ), the sys¬ 
tem evolves instead towards the momenta Pko. The spa¬ 
tial distribution appears to be identical to the uniform 
pumping case, but the phase changes by tt for adjacent 
maxima. This may be compared to the density and phase 
profiles in the ring geometry in Eig. which shows the 
same variation. 

This above shows that, a spatially dependent pump 
may be used to target particular momentum modes with 
the same periodicity. We note that without an initial seed 
momentum state, the pumping technique cannot gener¬ 
ate the desired momentum modes. The reason is that 
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FIG. 8: Schematic for sensitivity enhancement effect for pe¬ 
riodic potentials, (a) The real part of the wavefunction in 
a Bloch band (solid line) and the periodic potential (dashed 
line), (b) The resulting interference due to two waves of op¬ 
posing momentum. Due to the spatial modulation of the 
Bloch wavefunction, the resolution of the interference part- 
tern is enhanced. 
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FIG. 7: Metastable condensation dynamics, (a) The ini¬ 
tial state with equal populations of the three momenta k — 
0, /co, —k{). Time evolved steady state distributions for (b) uni¬ 
form pumping Eq. (13) and (c) periodic pumping Eq. (14). 
Parameters used are 7 = 77 = Vb = = 1, To = 2, /co = 27r7T0, 

^fe=o = ^k=kQ — ^k=-ko = 1/a/ 3- Timescales are in units of 
21710^Ih where a is the lengthscale unit. 


unlike the coherent excitation technique of Sec. II A[ the 
pumping profile P{x,y) does not contain any phase in¬ 
formation, hence cannot impose any particular phase on 
the condensate. Thus some population in the desired 
momentum mode must be there initially, but from there 
the pumping can reinforce this population, and diminish 
the other contributions. In practice, below the condensa¬ 
tion threshold many momentum modes are occupied and 
hence as only crosses the critical pump density the system 
should naturally gravitate towards the target momentum 
states. 


B. Periodic potentials 

Another strategy that can be used to enhance the sen¬ 
sitivity is via the use of periodic potentials on the un¬ 
derlying DBR lattice. The origin of this effect is the 
increased effective mass of the polaritons in the periodic 
lattice. In general, when a periodic potential is applied 
to a system, the energy-momentum dispersion breaks up 
into bands. Within each band, the stronger the potential 
that is applied, the flatter the dispersion becomes. This 
originates from the fact that the kinetic energy is sup¬ 
pressed due to the reduced tunneling between adjacent 
sites. Flatter bands result in a larger effective mass of 


the particles. As the de Broglie wavelength of the parti¬ 
cles is A = h/mv^ where v is the velocity of the particle, a 
large mass gives a shorter wavelength, and thus increased 
sensitivities. 

A more direct physical argument which shows the in¬ 
creased sensitivity is shown in Fig. Due to Bloch’s 
theorem, any eigenstate in a periodic potential may be 
written 


^Pkix) = e^’^'^Ukix) (16) 

where Uk{x) is a periodic function with the same periodic¬ 
ity as the crystal and k is the wavenumber. The envelope 
function Uk{x) can be decomposed into Wannier func¬ 
tions, which for strong potentials can be approximated 
by the localized potential around one of the potential 
minima in the periodic function. Thus in a periodic po¬ 
tential with counter-propagating modes the interference 
pattern is 

^ikox _j_ ^-ikox 

I-- Uka{x)f = cos^ koxul^ix). (17) 

As illustrated in Fig. the periodic potential thus gives 
a finer periodicity than the ±ko modes. By detecting 
shifts in the finer fringes this can give rise to higher 
sensitivities for the gyroscope. The mechanism for the 
increased sensitivity originates from the introduction of 
high momentum components in the periodic potential 
Uk{x). This is consistent with the more straightforward 
strategy of increasing /cq, except that the high momen¬ 
tum modes are naturally present in a periodic lattice. For 
the standard approach one must inject the system with 
a momentum T/cq. 

To check the stability of vortex-antivortex superposi¬ 
tions in periodic lattice potentials, we again perform a 
time evolution under Eq. <§• Specifically, we choose a 
Kagome or a basketweave lattice [59]. A Kagome lattice 
consists of interlaced triangles and exhibits a higher de¬ 
gree of frustration when compared with other 2D lattices 
(e.g., square and triangular lattices) [60]. A Kagome lat¬ 
tice is particularly interesting as in the tight-binding limit 
with only nearest-neighbor tunneling the 2D Kagome 
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lattice geometry possesses a completely flat band. In 
such flat bands, the polariton condensate order param¬ 
eters become tightly localized at the potential dips [6T] . 
Experimentally, several schemes have been proposed to 
create flat bands in the 2D Kagome lattice, utilizing 
metallo-photonic waveguides [62] , photonic crystal struc¬ 
tures [63] , or depositing a thin metal film [ST] |58l EH ES] • 
Flat bands can also be engineered in a frustrated lat¬ 
tice of micro-pillar optical cavities [66] . The lowest order 
Fourier decomposition of a Kagome lattice potential Vext 
can be described using just three components, written 

Ve^t = Vo , (18) 

where /i(x) = e*^opa^/(i+4p/3) cos(/copx/(l + 4p/3)), bi = 
(l/(l + 4p/3),0), b2 = (-l/(2 + 8p/3),-V3/2), bs = 
( — 1/(2 + 8p/3), V^/2), Vb > 0 is the lattice intensity, ko 
is a scalar quantity proportional to the lattice constant, 
and p ^ 3/2 [63(68]. We implement the above potential 
in the dGP equation Eq. (§, and evolve forward in time 
starting from a resonantly induced vortex-antivortex su¬ 
perposition state given by Eq. 

In Fig. [^we show the spatiotemporal evolution of the 
vortex superposition state in a Kagome lattice poten¬ 
tial. Upon examination of Fig. [^ we can see that in the 
presence of a Kagome lattice, our vortex superposition 
state is coherent and stable as it is generated in the dis¬ 
tributed Bragg reflector (DBR) microcavity outlined in 
Fig. [T] Fig. 1^ also illustrates the phase of the two- 
dimensional vortex superposition non-equilibrium con¬ 
densate in a Kagome lattice potential, which shows the 
characteristic phase difference between the lobes of tt. It 
should be pointed out that the flat band in the Kagome 
lattice corresponds to the third lowest energy band. In 
our simulations the polaritons predominantly occupy the 
first band, hence the only effect of the bands is to pe¬ 
riodically modulate the interference. In order to target 
the flat band, techniques such as metastable condensa¬ 
tion are required to create the condensate in these high 
momentum states. 


IV. SAGNAC PHASE 

A. Comparison to optical Sagnac interferometers 

Optical Mach-Zehnder interferometers [69] are typi¬ 
cally capable of achieving an excellent rotational mea¬ 
surement sensitivity, thanks to the large area that the in¬ 
terferometric loop circumscribes [701 m]. In this respect, 
atomic matter wave interferometers typically have small 
loop areas and short-time rotational sensitivities [72l[73]. 
However, matter wave interferometers have the advan¬ 
tage that their de Broglie wavelength is much shorter 
than the optical counterparts. For atomic BECs, this 
makes the rotational measurement sensitivity of atomic 
matter wave interferometers per unit area exceed that 



FIG. 9: Time evolution of a polariton BEC initially seeded 
with a vortex-antivortex superposition for a Kagome lattice 
potential given by Eq. (18). Timeframes for the (a) initial 
state t = 0; (b) intermediate state t = 5; and (c) steady-state 
t = 10 are shown. Timescales are in units of h/meV. The 
following parameters were used in the simulations: Tig — 0.05, 
Ti^ = 1.0, and hg = 0.1, where energies are in meV. 


of optical ones by the ratio mc^/Tiuj ^ 10^^ [36]. Typ¬ 
ically polariton masses are extremely insignificant (i.e., 
rueff ^ 10“^me), where rUe is the bare electron mass. 
This almost entirely negates the advantage of polari¬ 
tons to rriQf^c^/Tiuj ^ 1, which is the same level as pho¬ 
tons. This is not surprising as polaritons themselves are 
part-matter, part-photon quasiparticles, hence much of 
their characteristics are inherited from light. In terms of 
sensitivity, polaritons would appear to be ill-suited for 
such gyroscopic applications. We argue below that for 
counter-rotating vortex superposition based interferome¬ 
ters this argument does not directly apply, as the Sagnac 
phase is independent of the de Broglie wavelength. 

Let US first review the state-of-the-art optical Sagnac 
interferometers. Typically, in an optical interferometer, 
light of wavelength A propagates in a closed loop. For 
each completed propagation in the closed loop configura¬ 
tion, such as in a fiber optic gyroscope, the Sagnac phase 
is given by 




SttAiood^ 


( 19 ) 


Ac 




















10 


where the A is the wavelength of the light, Q is the angu¬ 
lar (rotational) velocity, c is the speed of light, and Aioop 
is the area of the Mach-Zender interferometer. In a ring 
laser gyroscope, the quantity that is measured is the beat 
frequency of two counterpropagating lasers. In terms of 
the phase difference, this is now a time-evolving quantity 
of the form 




Sir A\QQp^^t 
Xp 


( 20 ) 


where p is the perimeter of the interferometric loop. The 
fiber optic gyroscope and the ring laser gyroscope phases 
differ in operation in that the latter is time dependent, 
while the former is not. This means that as the sys¬ 
tem rotates, the fiber optic gyroscope produces a shift 
in the fringes only when there is a rotation 11, while the 
ring laser gyroscope detects the total time-varying phase 
of the rotation. Therefore the longer one waits between 
comparisons of the interference fringes, the larger the an¬ 
gle difference, and therefore the sensitivity. 

In both types of optical gyroscopes, the phase is pro¬ 
portional to the area of the interferometric loop. This 
is quite different to the situation in a BEC where the 
phase around a vortex is fixed purely by the topological 
winding number. Thus independent of the path, a phase 
of 27tI is picked up around a vortex. For a fixed /, the 
path can in principle enclose an arbitrary area and yet 
will not pick up any additional phase. It is also indepen¬ 
dent of the de Broglie wavelength A for the same reasons. 
Therefore, for a counter-rotating vortex-antivortex super¬ 
position, the Sagnac phase is independent of the area of 
the BEC and the mass of the particles (i.e., polaritons) 
involved. As the Sagnac phase directly originates from 
the interference of two counter-rotating angular momen¬ 
tum, it follows that the the Sagnac phase can then be 
written as 


0o(t) = 2/Ot, 


( 21 ) 


and can be determined from the O of the vortex super¬ 
position state in a polariton BEC. 

Is it possible to regain the geometry and wavelength 
dependence of Eq. (19) in a BEC to further improve 


the sensitivity? For the ring geometry considered in Sec. 

the area and wavelength dependence is reinstated, 
as we show below. For the BEC case the Sagnac phase 
for polaritons can be written m 




( 22 ) 


where N{t) is the number of times a polariton revolves 
around the loop in a time t. In a ring geometry of radius 
r as seen in IMl US, we have for counter-propagating 
momenta Eko 


N{t) = 


thko 

2iTrm 


This therefore gives 


= 2korQt. 


(23) 


(24) 


We note here that each polariton in the BEC cloud has 
the same magnitude of the angular momentum, which 
contributes to the signal-to-noise ratio as will be dis¬ 
cussed below. The relation now has a boost in the sensi¬ 
tivity with the radius of the ring, which arises from the 
fact that the velocity of a wave with momentum ko has 
a value v = hkojm. This is in contrast to a vortex which 
decays with radius v = hl/mr. 

In the ring geometry we regain the enhancement due to 
the de Broglie wavelength in Eq. ( [^ , as for heavy parti¬ 
cles, the typical wavenumber /cq = 27r/A is much higher, 
which also results in enhanced sensitivities. In this sense, 
the use of periodic potentials as discussed in the previ¬ 
ous section can aid the sensitivity, as this produces ef¬ 
fectively larger effective masses, and hence larger typical 
ko. We note that on setting k^ = 27r/A, Aioop = 7rr^, 
and p = 27rr, Eq. (24) reduces to Eq. (20), so that BEC 
based gyroscopes can be considered to be equivalent to 
the ring laser gyroscopes up to geometrical factors and 
the wavelength of the particles being used. 


B. Detection of Sagnac phase 

When the laboratory frame rotates in time, the two 
counter-rotating matter waves pick up different phases, 
and the transverse density profile at ^ = 0 is given as 

(1 + cos[2/(0 + Qt)]) |'0(x, t)|^ , (25) 

or for the case of the ring geometry, 

(1 + cos[2/cor((/) + fit)]) |'0(x,t)|^ , (26) 

signifying a rotation of the interference pattern as the 
Sagnac phase accumulates over time. Following the ro¬ 
tation of the lobes in the interference pattern, one may 
determine the angular velocity fl from the density pro¬ 
file of the polariton BEC. The condensate density will 
show an interference pattern determined by the phase 
difference between the amplitudes and charges of the two 
vortex components, as seen in Figs. i-H This allows 
the polariton BEC to be used for the measurement of 
changes in the Sagnac phase caused by the rotation of 
the laboratory frame of reference (i.e., the interferometer 
itself). This can be performed using existing experimen¬ 
tal techniques where a CCD camera is used to image the 
photoluminescence from the polariton BEC. 

Determining the Sagnac phase using our proposed po¬ 
lariton Sagnac interferometer will depend on the vortex 
charge I of the polariton BEC superposition state. In 
practice, spiral phase plates (SPPs) generally work best 
for low integer values of 1. However, in principle, I can 
take on any integer value, although low integer values 
of I generally give a more stable vortex superposition. 
The proposed polariton Sagnac interferometer is based on 
imaging the standing wave of the polariton BEC, which is 
comparably much smaller than the area occupied by com¬ 
parable optical technologies (e.g., ring-laser gyroscope). 
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Thus the polariton ring geometry has the potential for 
higher sensitivities. 

We are now in a position to compare the rotational 
phases per particle of the interference patterns (pQ. Due 
to the similar nature of the operation, it is simplest to 
compare the polariton gyroscopes to the ring laser gy¬ 
roscope, which both have an increasing Sagnac phase in 
the time domain. For the ring laser gyroscope, taking 
^loop ^ 1 m^ and p ^ 1 m, and typical optical wave¬ 
lengths, one obtains the proportionality of Qt in Eq. ( [^ 
to be SirAioop/ Xp ^ 10^. For ultracold atom implementa¬ 
tions, ^loop ^ 1 mm^, p 1 cm, and velocities v ^ Im/s 
we have SnAioo^/Xp ^ 10^ [76]. For the polariton gyro¬ 
scope with ring radii of r ^ 100/im and ko ^ 10pm~^ 
this gives the coefficient of the rotational phase in 
Eq. (24) as k^r ^ 10^. Thus in terms of the best achiev¬ 
able sensitivities per particle, polariton gyroscopes ap¬ 
pear to have a large handicap compared to other meth¬ 
ods due to their long wavelengths and relatively small 
interferometric loop areas. However, we shall see in the 
next section that some of this deficit is made up by the 
superior signal-to-noise of polaritons. 


C. Signal-to-noise ratio and stability 

We now discuss the signal-to-noise ratio (SNR) and the 
sensitivity achievable in such a polariton Sagnac interfer¬ 
ometer. The SNR is given by the ratio of the rotational 
phase shift to the shot-noise that depends on the num¬ 
ber of photons N arriving at the CCD from the polariton 
condensate per second [ZlITZl 

SNR=cPn{t)/VW- (27) 

The final SNR has a maximum at a certain rate of pho¬ 
ton detection, and the sensitivity f^min is calculated by 
setting SNR= 1. Thus, for the superposition in the ge¬ 
ometries as seen in Figs. i-E 

= l/i2ltVN), (28) 

and for the vortex-antivortex superposition in the ring 
geometry of Fig. [^we have 

= l/i2kortVN). (29) 

Typical polariton densities are in the region of 10^ cm“^, 
and a polariton spot size of 100 gives a polariton 
number in the region of 10^. The lifetime of the po¬ 
laritons is in the region of ^ 10 ps, hence we can esti¬ 
mate that the number of signal photons is N ^ 10^^ s“^ 
for polaritons. This gives our estimate for the polariton 
vortex-antivortex superposition in the ring geometry 

f^min ~ 10”^° rad s“^ (30) 

This is comparable with state-of-the-art atomic beam 
gyroscopes which are at the 10“^^ rad s~^ level 

[721176]. 


With regard to the stability of our proposed polariton 
Sagnac interferometer, it bypasses the beam drift that 
is a common cause of instability in the metrology that 
employs atomic or optical beams, and no complicated 
atomic beam configurations are required. The technical 
issues that optical Sagnac interferometers face, such as 
the short-time rotational sensitivity problem [78|, and 
the change in phase velocity of light m, should not 
be present for polaritons by construction. Moreover, 
thermal expansion over atomic beam trajectories is not 
an issue for a polariton BEC, as atomic beams are not 
needed. Interferometers which use atomic beams usually 
have beam drifts on the order of p deg /h after hours of 
running time [80]. The simplicity of our scheme (i.e., no 
large optical loops in a solid state device), would sug¬ 
gest that our scheme should not suffer as seriously from 
stability issues. 

We note that for purposes of seismometry, since the 
vector of angular velocity ft is proportional to the vector 
of ground velocity v 


ft = X V. (31) 

In this sense a quantum polariton seismometer can detect 
the ground velocity from the Sagnac phase by detecting 
rotations that originate from unaccounted rotations of 
the ground. 


V. SUMMARY AND CONCLUSION 

In this study, we have introduced two schemes for in¬ 
ducing counter-rotating vortex superposition states in 
a polariton BEC for gyroscopy. The first is via direct 
injection of orbital angular momentum (0AM) of light 
through a distributed Bragg reflector (DBR) microcavity. 
Once the vortex superposition state is induced, it follows 
a free evolution under standard pump-loss dynamics of 
the polariton BEC. The second is via metastable conden¬ 
sation in periodic lattices, where particular momentum 
modes are reinforced by using pumping profiles with the 
target periodicity. By directly numerically simulating the 
dynamics, we found that the vortex-antivortex superpo¬ 
sition state is stable under a variety of different condi¬ 
tions; in smooth, disordered, and periodic potentials. In 
all cases we found that the vortex superposition state is 
stable and is long-lived. 

We investigated the possibility that the superposition 
of counter-rotating currents in the polariton BEC can be 
used as a method for determining the Sagnac phase - and 
thus the angular velocity of the interferometer - effec¬ 
tively constituting a “quantum gyroscope.” The vortex- 
antivortex superposition has the same Sagnac phase de¬ 
pendence as atomic BECs, as the phase relation around 
a vortex is topologically fixed, regardless of the particle 
mass or the area of the BEC. Optical Sagnac interfer¬ 
ometers have large sensitivities by increasing the area of 
the interference loop, and atomic Sagnac interferometers 
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have an advantage due to their small de Broglie wave¬ 
length. In this respect, polariton gyroscopes have rela¬ 
tively small loop areas and long wavelengths. This can 
be enhanced using ring geometries or periodic potentials. 
In comparison with current technology it would appear 
that on a per polariton basis atomic and optical interfer¬ 
ometers have a larger sensitivity. However, the polaritons 
do possess the great advantage that their signal-to-noise 
ratio should be considerably larger than that of atomic 
BECs. In an atomic BEC the optical measurement has 
the limitation that it should not disturb the system [81] , 
hence there are restrictions on how bright the probe light 
should be. In contrast, for polaritons they are themselves 
part light, hence illumination by laser light is intrinsic to 
the condensation process. Thus while the phase depen¬ 
dence of the polariton Sagnac interferometer starts with 
several orders of magnitude deficit, the final estimated 
sensitivity is competitive with existing methods. 

Polariton BEC based gyroscopes have the additional 
practical advantage to atomic BECs that they should 
lead to compact and room temperature compatible de¬ 
vices, as it is based on semiconductor technology. This 
could lead to a device with the desirable characteristics 
of compact operation with very high sensitivity and sta¬ 
bility. Eor the metastable condensation technique, the 
ability of electrical injection of the condensate would re¬ 
move the need for optics entirely, simplifying the device 
further. Due to the limitations in the size of the polari¬ 
ton condensate, they are naturally suited towards small 
volume applications. The effective mass of the polaritons 
can be tuned by engineering periodic structures, and con¬ 
ventional methods such as changing the exciton fraction 
of the polaritons, which give additional boosts to sensi¬ 
tivity. 
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Appendix A: Orbital Angular Momentum states of 

light 


0AM superposition states of light can be generated 
by passing a TEMqo beam mode through a spiral phase 
plate (SPP), followed by a Mach-Zehnder interferometer 
with a Dove prism. The SPP converts the TEMqo to a 
Laguerre-Gaussian (i.e., LGi^p) beam mode of quantized 
0AM with quantum numbers I and p, and the Dove prism 
changes the winding of the quantized 0AM from Ih to 
—Ih. In addition to generation using spiral phase plates 
[82] , cylindrical lens mode converters [83] , and computer 
generated holograms [84] are alternative methods to gen¬ 
erate integer quantized 0AM. The LGi p beam modes 
form a complete orthonormal basis set of solutions for 
paraxial light beams commonly found in lasers. These 
modes can be expressed as [85] 


ul,p{r, 4>, z) 

X W I 


/ 2p! 

1_ 

\i\ 

_p2 

y ww'^{z){\i\+p)\ 

w{z) 

exp 



2r2 \ r 
w^(z) J 



(p{z) 


(AI) 


where the radius of the beam squared is w‘^{z) = 2(z^ + 
b‘^)/kb^ the radius of curvature is R{z) = {z^ + 
the Gouy phase ip{z) = (2p + |/| + I) tan“^(2;/6), b is 
the Rayleigh range, k is the wave number, and Lp^{r) 
describe the Laguerre polynomials. 


4''w = E^- 


i-iy 


{i+p)\ 


m=0 


{p — m)!(|/| + m)\m\ 


(A2) 


In Eqs. ( |AI| ) and ( |A2| ), p is the number of nonaxial radial 
nodes, and the index I is known as the winding number. 
The winding number describes the helical structure of 
the wavefront and the number of times the phase jumps 
occur around the beam path in the azimuthal direction. 
When / = p = 0, the LGi^p beam is identical to the 
fundamental Gaussian beam, or TEMqq. 

The higher order modes from the fundamental mode 
may be calculated using the operator algebra formalism 
[86] . Analogous to the two-dimensional quantum har¬ 
monic oscillator ISZIIHH] , at z = 0, 


i.(0) = 


\/2bk 


{kr cos((/)) + bcos{(l))dr 


and 


br sin((/))a0), 
(A3) 


Ay{0) = _ (fer sin(0) + 6sin(0)5y. + cos(0)50). 

V2bk 

(A4) 

These operators can be evolved in the 2;-direction by a 
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propagator U{z) for = 5^9^ according to 


The operators in Eq. (A8) obey the commutation rules 


-^P^z 

2k 


U{z) = exp 
= exp 

Next, we can rewrite Eqs. (|A3|) and (|A4[) as 


-^id^r+r^dl)z 


A,{z) = U{z)A,{Q)U\z), 
Ay{z) = U{z)Ay{0)U\z), 


(A5) 


A±{z),Al^{z) 


1, A±{z),A]^{z) 


= 0 . 


(A9) 


(A6) 


Higher-order LGi^p beam modes can be obtained from 
operating on the TEMqo according to [87] 


(A7) 


where 


1 

71 

1 


2\/Ui 


A^{z) iAy{z) 


[kre^^‘i‘ + + 2ze^^^){dr + r^^)] . 


(A8) 


= {^\ [AwJ'tEMooW, 


(AlO) 

where m = (/+p)/2 and n = {l—p)/2. We note that 0AM 
states ufp{z) and u^p{z) differ only in the direction of 
the phase winding, that is, clockwise or counterclockwise, 
respectively. 
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